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ABSTRACT 

We report on a correlation between virial mass M and spin parameter A for dark matter halos forming at redshifts 
z ^ 10. We find that the spin parameter decreases with increasing halo mass. Interestingly, our analysis indicates 
that halos forming at later times do not exhibit such a strong correlation, in agreement with the findings of previous 
studies. We briefly discuss the implications of this correlation for galaxy formation at high redshifts and the galaxy 
population we observe today. 

Subject headings: galaxies: formation — cosmology: theory — cosmology: early Universe — methods: 
numerical 



L INTRODUCTION 

The physical mechanism by which galaxies acquire their an- 
gular momentum is an important problem that has been the sub- 
ject of investigation for nearly sixty years (Hoyle 1949). This 
reflects the fundamental role played by angular momentum of 
galactic material in defining the size and shapes of galaxies (e.g. 
Fall & Efstathiou 1981). Yet despite its physical significance, 
a precise and accurate understanding of the origin of galactic 
angular momentum remains one of the missing pieces in the 
galaxy formation puzzle. 

A fundamental assumption in current galaxy formation mod- 
els is that galaxies form when gas cools and condenses within 
the potential wells of dark matter halos (White & Rees 1978). 
Consequently it is probable that the angular momentum of the 
galaxy will be linked to the angular momentum of its dark mat- 
ter halo (e.g. Fall & Efstathiou 1980; Mo, Mao & White 1998; 
Zavala, Okamoto & Frenk 2007). Within the context of hier- 
archical structure formation models, the angular momentum 
growth of a dark matter proto-halo is driven by gravitational 
tidal torquing during the early stages (i.e. the linear regime) of 
its assembly. This "Tidal Torque Theory" has been explored in 
detail; it is a well-developed analytic theory (e.g. Peebles 1969, 
Doroshkevich 1979, White 1984) and its predictions are in 
good agreement with the results of cosmological A^-body sim- 
ulations (e.g. Barnes & Efstathiou 1987; Warren et al. 1992; 
Sugerman, Summers & Kamionkowski 2000; Porciani, Dekel 
& Hoffman 2002). However, once the proto-halo has passed 
through maximum expansion and the collapse has become 
non-linear, tidal torquing no longer provides an adequate de- 
scription of the evolution of the angular momentum (White 
1984), which tends to decrease with time. During this phase it 
is likely that merger and accretion events play an increasingly 
important role in determining both the magnitude and direction 
of the angular momentum of a galaxy (e.g. Bailin & Steinmetz 
2005). Indeed, a number of studies have argued that mergers 
and accretion events are the primary determinants of the an- 
gular momenta of galaxies at the present day (Gardner 2001; 
MaUer, Dekel & Somerville 2002; Vitvitska et al. 2002). 

It is common practice to quantify the angular momentum of 
a dark matter halo by the dimensionless "classical" spin param- 
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eter (Peebles 1969), 



A = 



(1) 



where J is the magnitude of the angular momentum of material 
within the virial radius, M is the virial mass, and E is the total 
energy of the system. It has been shown that halos that have 
suffered a recent major merger will tend to have a higher spin 
parameter A than the average (e.g. Hetznecker & Burkert 2006; 
Power, Knebe & Knollmann 2008). Therefore one could argue 
that within the framework of hierarchical structure formation 
that higher mass halos should have larger spin parameters on 
average than less massive systems because they have assem- 
bled a larger fraction of their mass (by merging) more recently. 

However, if we consider only halos in virial equilibrium, 
should we expect to see a correlation between halo mass and 
spin? One might naively expect that more massive systems will 
have had their maximum expansion more recently and so these 
systems will have been tidally torqued for longer than systems 
that had their maximum expansion at earlier times. This sug- 
gests that spin should increase with timing of maximum expan- 
sion and therefore halo mass. However, one finds at best a weak 
correlation between mass and spin for equilibrium halos at z=Q 
(e.g. Cole & Lacey 1996; Maccio et al. 2007; Bett et al. 2007, 
hereafter B07), and the correlation is for spin to decrease with 
increasing halo mass, contrary to our naive expectation. 

In this paper, we report on a (weak) correlation between spin 
and mass for equilibrium halos at redshift z=10. The trend is 
for higher-mass halos to have smaller spins, and is qualitatively 
similar to the one reported by B07 for the halo population at 
z=0. We present the main evidence in support of this correla- 
tion in Section 4 and we consider its implications for galaxy 
formation in Section 6. 

2. THE SIMULATIONS 

For the simulations presented in this paper we have adopted 
the cosmology as given by Spergel et al. (2003) (ilo = 0.3, 
flA = 0.7, (Tg = 0.9, and Hq = 70km/sec/Mpc). Each run em- 
ployed = 256^ particles and differed in simulation box-size 
Lbox, which leads to the particle mass differing between runs 
- nip = p„ni^o(L\,ox / N)^ , where pcrit = 3Hq/SttG. This allows 
us to probe a range of halo masses at redshift z= 1 0. The primary 
parameters of these simulations are summarized in Table 1 . 
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Halos in all runs have been identified using the MPI paral- 
lelized version of the AHF halo finder-' (AMIGA's-Halo-Finder), 
which is based on the MHF halo finder of Gill, Knebe & Gibson 
(2004). For each halo we compute the virial radius R, defined 
as the radius at which the mean interior density is Avir times the 
background density of the Universe at that redshift. This leads 
to the following definition for the virial mass M: 



M: 



47r 3 

AvirSipcrit^ 



(2) 



Note that Avir is a function of redshift and amounts to Av 



210 at redshift z=10 , Av 



230 at z=l, and the "usual" 



Avir « 340 at z=0 (cf. Gross 1997). Table 1 sunamarises the to- 
tal number of halos (A'haios) recovered by AHF, while Zfinai gives 
the redshift that the simulation has been evolved to. 

We add that we ran five reahsations of B20 to redshift z=10 
in order to have a statistically significant sample of halos in that 
particular model. However, we also note that the fitting param- 
eters presented in the following Section are robust in the sense 
that they do not depend on whether we stack the halos from 
those five runs or use them individually. 

3. THE HALO SAMPLE 

We show in a companion paper (Power, Knebe & Knollmann 
2008) that a substantial fraction of the halo population at high 
redshift is not in virial equilibrium. Because we wish to exam- 
ine the spin distribution of equilibrium halos, it is important to 
account for unrelaxed systems when investigating correlations 
between spin and halo mass. For example, it has been shown 
that the spin can increase sharply in the aftermath of mergers 
with mass ratios as modest as 5:1 (e.g. Hetznecker & Burkert 
2006), and that the degree to which a halo is in dynamical relax- 
ation is as important as recent merging history in its influence 
on spin (D'Onghia & Navarro 2007). To ensure that halos in 
our sample are in virial equilibrium, we compute the virial ratio 
for each halo, which we define as 



Q 



2T + S„ 



U 



+ 1 



(3) 



Here T represents the kinetic energy, U the potential energy, 
and Sp the surface pressure of a given halo of mass M. By in- 
cluding Sp, we can account for the effect of infalling material 
on the dynamical state of the halo. Each of these quantities 
are evaluated using all gravitationally bound particles, and we 
adopt the formula of Shaw et al. (2006) for the surface pressure 
term Sp (cf. equations. (4)-(6) in their study). 

In Figure 1 we show that the relation between halo mass and 
Q can vary with mass. This is apparent at redshift z=l, where 
we find a trend for more massive halos to be less viriaUsed. In 
contrast, high redshift halos are less virialised on average (as 
indicated by the increased average {Q) « —0.3), but we find no 
apparent trend with mass. 

Why is there a mass dependence at z=l but not at z=10? 
There are two factors. The first is that high redshift halos "see" 
an effective slope of the initial power spectrum of ngff ~ -3, and 
so the time at which a particular mass scale starts to collapse 
is relatively insensitive to mass. Therefore we do not expect 
to find a strong correlation between virial state Q and mass. 
We have checked this halo populations drawn from the simula- 
tions of scale-free cosmologies of Knollmann, Power & Knebe 



(2008) and our interpretation is consistent with the correlations 
we find in these runs. The second is that the typical collapsing 
mass M* at z=10 is small - of order lO^/i'^M© - and because 
we resolve mass scales that have collapsed more or less simul- 
taneously, we see a population that has yet to relax. At z=l, the 
typical collapsing mass is much larger - of order lO"/z~'M0 - 
and so we resolve a population of halos whose mass accretion 
histories are more diverse. The most massive systems tend to 
be ones that have formed most recently, and are therefore the 
least dynamically relaxed. 

We have used the following relation between Q and M to 
classify dynamically relaxed and unrelaxed systems: 



2allowed 
Gallowed 



a 
oc 



^-0.015 

const. 



, for z = 1 
, for z = 10 



(4) 



We allow the Q values of halos in our sample to deviate from 
these scaling relations by not more than 



allowed ' 



■Glim<G<G. 



illowed 



lim 



(5) 



with Glim = 0.15 (indicated by the dashed lines in Fig. 1). Fur- 
thermore, we consider only halos that contain at least A^min = 
600 particles within their virial radius to ensure that we are not 
influenced by particle discreteness. Interestingly, when com- 
puting spin, the tighest restriction on particle number comes 
not from the calculation of angular momentum but from the 
calculation of the potential energy. By comparing analytic solu- 
tions with Monte Carlo realisations of Navarro, Frenk & White 
(1997) haloes, we find that at least 600 particles are required if 
the energy is to be computed to better than 10%. 

4. THE SPIN-MASS CORRELATION 

Calculating the total energy £ of a halo is computationally 
expensive, and so computing A using Eq. (1) is also expensive. 
This prompted Bullock et al. (2001) to introduce a modified 
spin parameter 



A' = 



J 



(6) 



V2MVR ' 

where V = \/GM/R measured the circular velocity at the virial 
radius R and J represents the absolute value of the angular mo- 
mentum. We follow Bullock et al. and compute spin using 

equation (6). 

In Figure 2, we investigate the correlation between halo spin 
A' and mass M. We show only halos that fulfill our selec- 
tion criteria (individual dots) and bin the data in five mass bins 
equally spaced in log-space between Mmi,, and M^ax of the con- 
sidered halos at the respective redshift. The values plotted as 
histograms thereby represent the weighted mean of all spin pa- 
rameters in the respective mass range where the weight is in- 
versely proportional to the error estimate 



J 



0.2 



(7) 



for the spin parameter of a halo consisting of A' particles as de- 
rived in Bullock et al. (2001) (cf. equation (7) in that study). 
The error bars indicate the standard deviation of the spin pa- 
rameter values in the bin from the weighted mean."* 

The best fitting power-laws to these histograms reveal that 
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A' cx M« 



with 



a = -0.002 ±0.149 ,forz=l 
a = -0.059±0.171 ,forz=10. 



(8) 



(9) 



This indicates that there is a weak correlation at high red- 
shifts for spin to decrease with increasing mass, albeit stronger 
than the one at z=l. We compute Spearman rank correlation 
coefficients at z=10 (1) and find /?, = -0.137(0.06). 

As an alternative approach, we fit a lognormal function to 
each of our halo samples at z=10 and z=l, 



P(A') = 



1 



A'v27r (To 



exp I 



2al 



(10) 



The resulting curves are presented in Fig. 3 whereas the best-fit 
parameters, median values for A^g^ = median(A') and median 
halo masses M^ed are given in Table 2. Inspection of the best- 
fit parameters confirm that the median spin declines as we move 
from less massive to more massive objects at high redshift. 

5. STABILITY OF RESULTS 

Because of the weak nature of the measured correlation it is 
vitally important to check its credibihty by performing a statis- 
tical analysis. To this extent we investigate the sensitivity of the 
logarithmic slope a with respect to a number of parameters that 
enter into its determination, namely the number of bins A^wns 
used for the histograms; the virialisation criterion parametrized 
via 2iim; and the minimum number of particles AWn within a 
halo's virial radius. Note that we vary one parameter at a time, 
keeping the others at their fixed "standard" values. The results 
are presented in Table 3, Table 4, and 5. 

We find that bin number has practically no effect on the slope. 
Similarly we find that varying the viriahsation criterion has 
little effect on the slope of the relation between mass and spin, 
regardless of redshift (Table 4). In contrast, we find that the 
minimum number of particles within a halo's virial radius has 
a strong and systematic effect on the result at z=l - as A'min in- 
creases, we find that the logarithmic slope becomes shallower. 
This does not appear to be true at high redshift, although as we 
go to earlier times we find that the number of massive haloes 
becomes progressively smaller and our determination of a be- 
comes increasingly unreUable. 

These tests lead us to believe that our results are both stable 
and reliable, and our main result holds: the correlation between 
spin paramater A' and halo mass M is one order of magnitude 
larger at redshift z=10 than atz=l. 

As a further test of the credibility of the correlation we mea- 
sure between halo mass and spin, we use the criteria of three 
other studies to select our halo sample. These are: 

• Maccio et al. (2007) criteria: 

- N^,, = 25Q 

- Xos < 0.04 

- Prms < 0.4 

• Bett et al. (2007) criteria: 

- AU,=300 



- Glim = 0.5 

• Neto et al. (2007) criteria: 

- AUn = 600 

- Xoff < 0.07 

- /sub < 0.1 

Here A'min is again the minimum number of particles in a 
halo, Q\im the virial limit as defined in Eq. (5), JCoff measures the 
distance between the most bound particle and the centre of mass 
in units of the virial radius, prms is an indicator of how well the 
density profile of the halo can be fitted by a Navarro, Frenk & 
White (1997) profile (cf. Eq. (2) in Maccio et al. 2007), and 
/sub is the fraction of mass in subhalos. The resulting power 
law slopes a (cf. Eq. (8)) for z = 1 and z = 10 are presented in 
Table 6. Again we note that there appears to be a much stronger 
correlation between A' and M at higher redshift. While the re- 
lation is consistent with zero at z=l (as confirmed by Maccio 
et al. 2007) spin and mass are correlated at z=10. How this 
result relates to the Bett et al. (2007) result - who find a weak 
correlation at z=0 - will be discussed in the following Section. 

6. CONCLUSIONS AND DISCUSSION 

We have performed a careful investigation of the relation 
between virial mass and dimensionless spin parameter for dark 
matter halos forming at high redshifts z > 10 in a ACDM cos- 
mology. The result of our study, which is based on a series of 
cosmological A^-body simulations in which box size was var- 
ied while keeping particle number fixed, indicates that there is 
a weak correlation between mass and spin at z=10, such that 
the spin decreases with increasing mass. If there is a corre- 
lation at z=l, we argue that it is significantly weaker than the 
one we find at z=10; this is in qualitative agreement with the 
findings of previous studies that focused on lower redshifts 
(Maccio et al. 2007, Shaw et al. 2005, Lemson & Kauffmann 
1999). 

Interestingly, B07 find a weak correlation between me- 
dian spin and halo mass at z=0 in the Millennium Simulation 
(Springel et al. 2005), in the same spirit as the one presented 
here for z=10: lower mass halos tend to have higher spins. 
However, as we show, the correlation between halo mass and 
spin is weaker at z=l than at z=10, whereas the correlation re- 
ported in B07 for halos at z=0 is much stronger than the one 
we find at z=l. This is not what one would expect, and so 
it is important to try and understand the source of the differ- 
ence between our result at z=l and the B07 result at z=0. B07 
fitted a 3''''-order polynomial to the median spins of halos in 
the mass range 3 x 10" M/ih-^Mo) ^ 3 x 10"* at z=0. The 
form of this polynomial is extremely sensitive to the precise 
values of the best-fit parameters (Bett, private communication) 
and it is not straightforward to extrapolate its behaviour out- 
side of the given mass range and redshift. We derive our es- 
timates of the power-law exponents from the spin distribution 
with respect to halo mass at z= 1 . Our halos he in the mass range 
3x10^^ M/ih-^Mo) ^ 5 X lO^^. B07 base their median spins 
upon ^1.5 million halos with correspondingly small errors, 
and note that the weak nature of the trend of spin with mass 
makes it hard to detect. This suggests to us that the correlation 
between mass and spin at z=10 is remarkably strong rather than 
the correlation at z=l being too weak! 
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When studying correlations between halo mass and spin, 
great care must be taken in defining the halo sample. In partic- 
ular, we find that mass resolution (i.e. the number of particles 
with which a halo is resolved) and the degree of virialisation 
of a halo can have a significant effect on the strength of the 
correlation (at least at z=l, cf. Table 5). This - at least - is in 
good agreement with the findings of B07. 

We note that Power & Knebe (2006) demonstrated that the 
size of simulation box can lead to a suppression of angular mo- 
mentum in smaller boxes, due to the absence of longer wave- 
length perturbations in the initial conditions. This will lead 
to a bias in our estimate of A (approximately a ~ 10% effect) 
but we have verified that the spin distributions we obtain from 
a B20 run truncated on scales larger than the longest wave- 
length perturbation modelled in the Bl run produces results 
that are consistent. Indeed, we would expect the correlation 
to be strengthened if the B I spins were corrected for box size 
effects. 

It is interesting to speculate on the consequences of this cor- 
relation for galaxy formation at high redshifts and the galaxy 
population we observe today. In the standard picture of galaxy 
formation, gas cools on to dark matter halos and is shock heated 
to the virial temperature of the halo. The angular momentum of 
the gas and the dark matter should (initially) be similar because 
they are subject to the same tidal field. As the innermost dens- 
est parts of the gaseous halo cool, they will settle into a gaseous 



disk with a scale length determined by the specific angular mo- 
mentum of the gas, which we would expect to be related to the 
angular momentum of the halo (e.g. Zavala, Okamoto & Frenk 
2007). 

If more massive halos at high redshifts show a tendency to 
have smaller spin parameters, the gas disks will have lower spe- 
cific angular momenta and therefore will be more centrally con- 
centrated. If star formation rate correlates with surface density, 
then we might expect the star formation rate to be enhanced in 
more massive halos. Because massive halos tend to form prefer- 
entially in high density, highly clustered environments in which 
the merger rate also tends to be enhanced, then we might expect 
star formation to proceed more rapidly and at earlier times in 
these halos. Might this explain the effect of "downsizing" (e.g. 
Cowie et al. 1996), the successive shifting of star formation 
from high- to low-mass galaxies with decreasing redshift? We 
shall pursue this in a more quantitative manner in future work. 
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tiative", grant DP 0665574. The simulations presented in this 
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Astrophysics & Supercomputing, Swinburne University as well 
as the Sanssouci cluster at the Astrophysikalisches Institut Pots- 
dam. 
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Table 1 

Summary of the cosmological simulations and number of halos. 



ran 


Lbox [fr^ Mpc] 


ntp [h 1 Mq] 


Zfinal 


'^halos 


'^halos 


relaxed halos 


ArZ=10 

relaxed halos 


BOl 


1 


4.9 xlO^ 


10 




8780 




286 


B02 


2.5 


7.8 X 10"* 


10 




7991 




201 


BOS 


5 


6.2 X 10^ 


1 


16917 


6532 


832 


109 


BIO 


10 


4.9 X 10^ 


1 


18589 


4360 


949 


37 


B20 


20 


4.0 X 10^ 


1 


20514 


10947 


995 


27 


BlOO 


100 


4.9 X 10' 


1 


24696 




978 





Table 2 
Fitting parameters for PC A). 



ran 


K 


a 


\' 


M„,ed [h ^ Mq] 


ran 


K 


a 


'^med 


M„,ed [h ^ Mq] 








z=l 








z- 


= 10 














BOl 


0.042 


0.538 


0.041 


5.95e+06 












B02 


0.040 


0.539 


0.037 


8.19e+07 


B05 


0.033 


0.521 


0.031 


7.92e+08 


B05 


0.036 


0.516 


0.035 


5.63e+08 


BIO 


0.037 


0.527 


0.034 


6.15e+09 


BIO 


0.033 


0.540 


0.029 


4.14e+09 


B20 


0.038 


0.543 


0.036 


5.40e+10 


B20 


0.030 


0.251 


0.027 


3.33e+10 


BlOO 


0.037 


0.535 


0.035 


5.25e+12 













Table 3 

Variation of a with A^bins (Giim = 0.15, Nmia = 600). 







az=io±aa 


4 


-0.001 ±0.169 


-0.061 ± 0.191 


5 


-0.002 ±0.149 


-0.059 ±0.171 


6 


-0.005 ±0.137 


-0.053 ±0.155 


7 


-0.006 ±0.128 


-0.056 ±0.144 


8 


-0.003 ±0.120 


-0.069 ±0.132 



Table 4 

Variation of a with gum (A^bins = 5,AUn = 600). 



Glim 


a~i ± act 




WZ=1 
relaxed halos 


ArZ=10 

relaxed halos 


0.05 


-0.007 ±0.148 


-0.058 ±0.153 


2123 


264 


0.10 


-0.004 ±0.148 


-0.062 ±0.168 


3365 


497 


0.15 


-0.002 ±0.149 


-0.059 ±0.171 


3811 


660 


0.20 


-0.005 ±0.151 


-0.052 ±0.173 


3994 


739 


0.25 


-0.002 ± 0.153 


-0.052 ±0.173 


4165 


774 



Table 5 

Variation of o with A^mj,, (Qiim = O.lS.A^bins = 5). 



^min 


az=i ± (Ja 




relaxed halos 


relaxed halos 


100 


0.003 ± 0.129 


-0.037 ± 0.158 


19707 


5478 


200 


0.002 ±0.137 


-0.038 ±0.167 


10809 


2534 


300 


0.002 ±0.142 


-0.035 ±0.174 


7336 


1526 


600 


-0.002 ±0.149 


-0.059 ±0.171 


3811 


660 


1000 


-0.004 ±0.155 


-0.058 ±0.196 


2303 


343 


2000 


-0.011 ±0.164 


-0.048 ±0.186 


1154 


120 



Table 6 

Applying different virialisation criterion (A/bins = 5). 



criterion 


CKz=l 




relaxed halos 


relaxed halos 


Netoetal. (2007) 


0.000 ±0.149 


-0.041 ±0.178 


3486 


429 


Maccio et al. (2007) 


0.001 ± 0.146 


-0.040 ±0.184 


4275 


512 


Bett et al. (2007) 


0.003 ±0.146 


-0.040 ± 0.175 


8582 


1955 
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Fig. 1 . — Relation between virial parameter Q (cf. equation 3) and halo mass. The solid lines represent the adopted virialisation criteria as given by Eq. (4). Note 
that we already applied the mass cut of 600 particles per halo for this plot and hence the number of halos appearing does not agree with the number given in Table 1 . 
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Fig. 2. — CoiTelation of spin parameter A with mass M. The binned data (histograms) has been fitted to a power laws (dashed line, cf. Eq. (9)). 




Fig. 3. — Lognormal distributions of the spin parameter A'. 



